From NMR to AI: Designing a Novel Chemical Representation to Enhance Machine Learning Predictions of Physicochemical Properties.

A novel approach to the utilization of nuclear magnetic resonance (NMR) spectroscopy data in the prediction of logD through machine learning algorithms is shown. In the analysis, a data set of 754 chemical compounds, organized into 30 clusters, was evaluated using advanced machine learning models, such as Support Vector Regression (SVR), Gradient Boosting, and AdaBoost, and comprehensive validation and testing methods were employed, including 10-fold cross-validation, bootstrapping, and leave-one-out. The study revealed the superior performance of the Bucket Integration method for dimensionality reduction, consistently yielding the lowest root mean square error (RMSE) across all data sets and normalization schemes. The SVR prediction models demonstrated remarkable computational efficiency and low cost, with the best RMSE value reaching 0.66. Our best model outperformed existing tools like JChem Suite's logD Predictor (0.91) and CplogD (1.27), and a comparison with traditional molecular representations yielded a comparable RMSE (0.50), emphasizing the robustness of our NMR data integration. The widespread availability of NMR data in pharmaceutical and industrial research presents an untapped resource for predictive modeling, highlighting the need for accessible methodologies like ours that complement the analytical toolbox beyond conventional 2D approaches. Our approach, designed to leverage the rich spatial data from NMR spectroscopy, provides additional insights and enriches drug discovery and computational chemistry with a freely accessible tool.

Superiority of the Triple-Acting 5-HT6R/5-HT3R Antagonist and MAO-B Reversible Inhibitor PZ-1922 over 5-HT6R Antagonist Intepirdine in Alleviation of Cognitive Deficits in Rats.

The multifactorial origin and neurochemistry of Alzheimer's disease (AD) call for the development of multitarget treatment strategies. We report a first-in-class triple acting compound that targets serotonin type 6 and 3 receptors (5-HT-Rs) and monoamine oxidase type B (MAO-B) as an approach for treating AD. The key structural features required for MAO-B inhibition and 5-HT6R antagonism and interaction with 5-HT3R were determined using molecular dynamic simulations and cryo-electron microscopy, respectively. Bioavailable PZ-1922 reversed scopolamine-induced cognitive deficits in the novel object recognition test. Furthermore, it displayed superior pro-cognitive properties compared to intepirdine (a 5-HT6R antagonist) in the AD model, which involved intracerebroventricular injection of an oligomeric solution of amyloid-ß peptide (oAß) in the T-maze test in rats. PZ-1922, but not intepirdine, restored levels of biomarkers characteristic of the debilitating effects of oAß. These data support the potential of a multitarget approach involving the joint modulation of 5-HT6R/5-HT3R/MAO-B in AD.

Dual Piperidine-Based Histamine H3 and Sigma-1 Receptor Ligands in the Treatment of Nociceptive and Neuropathic Pain.

In search of new dual-acting histamine H3/sigma-1 receptor ligands, we designed a series of compounds structurally based on highly active in vivo ligands previously studied and described by our team. However, we kept in mind that within the previous series, a pair of closely related compounds, KSK67 and KSK68, differing only in the piperazine/piperidine moiety in the structural core showed a significantly different affinity at sigma-1 receptors (σ1Rs). Therefore, we first focused on an in-depth analysis of the protonation states of piperazine and piperidine derivatives in the studied compounds. In a series of 16 new ligands, mainly based on the piperidine core, we selected three lead structures (3, 7, and 12) for further biological evaluation. Compound 12 showed a broad spectrum of analgesic activity in both nociceptive and neuropathic pain models based on the novel molecular mechanism.

Tuning the Biological Activity of PI3Kδ Inhibitor by the Introduction of a Fluorine Atom Using the Computational Workflow.

As a member of the class I PI3K family, phosphoinositide 3-kinase δ (PI3Kδ) is an important signaling biomolecule that controls immune cell differentiation, proliferation, migration, and survival. It also represents a potential and promising therapeutic approach for the management of numerous inflammatory and autoimmune diseases. We designed and assessed the biological activity of new fluorinated analogues of CPL302415, taking into account the therapeutic potential of our selective PI3K inhibitor and fluorine introduction as one of the most frequently used modifications of a lead compound to further improve its biological activity. In this paper, we compare and evaluate the accuracy of our previously described and validated in silico workflow with that of the standard (rigid) molecular docking approach. The findings demonstrated that a properly fitted catalytic (binding) pocket for our chemical cores at the induced-fit docking (IFD) and molecular dynamics (MD) stages, along with QM-derived atomic charges, can be used for activity prediction to better distinguish between active and inactive molecules. Moreover, the standard approach seems to be insufficient to score the halogenated derivatives due to the fixed atomic charges, which do not consider the response and indictive effects caused by fluorine. The proposed computational workflow provides a computational tool for the rational design of novel halogenated drugs.

Impact of the Substitution Pattern at the Basic Center and Geometry of the Amine Fragment on 5-HT6 and D3R Affinity in the 1H-Pyrrolo[3,2-c]quinoline Series.

Salt bridge (SB, double-charge-assisted hydrogen bonds) formation is one of the strongest molecular non-covalent interactions in biological systems, including ligand-receptor complexes. In the case of G-protein-coupled receptors, such an interaction is formed by the conserved aspartic acid (D3.32) residue and the basic moiety of the aminergic ligand. This study aims to determine the influence of the substitution pattern at the basic nitrogen atom and the geometry of the amine moiety at position 4 of 1H-pyrrolo[3,2-c]quinoline on the quality of the salt bridge formed in the 5-HT6 receptor and D3 receptor. To reach this goal, we synthetized and biologically evaluated a new series of 1H-pyrrolo[3,2-c]quinoline derivatives modified with various amines. The selected compounds displayed a significantly higher 5-HT6R affinity and more potent 5-HT6R antagonist properties when compared with the previously identified compound PZ-1643, a dual-acting 5-HT6R/D3R antagonist; nevertheless, the proposed modifications did not improve the activity at D3R. As demonstrated by the in silico experiments, including molecular dynamics simulations, the applied structural modifications were highly beneficial for the formation and quality of the SB formation at the 5-HT6R binding site; however, they are unfavorable for such interactions at D3R.

Isomeric Activity Cliffs-A Case Study for Fluorine Substitution of Aminergic G Protein-Coupled Receptor Ligands.

Currently, G protein-coupled receptors (GPCRs) constitute a significant group of membrane-bound receptors representing more than 30% of therapeutic targets. Fluorine is commonly used in designing highly active biological compounds, as evidenced by the steadily increasing number of drugs by the Food and Drug Administration (FDA). Herein, we identified and analyzed 898 target-based F-containing isomeric analog sets for SAR analysis in the ChEMBL database-FiSAR sets active against 33 different aminergic GPCRs comprising a total of 2163 fluorinated (1201 unique) compounds. We found 30 FiSAR sets contain activity cliffs (ACs), defined as pairs of structurally similar compounds showing significant differences in affinity (≥50-fold change), where the change of fluorine position may lead up to a 1300-fold change in potency. The analysis of matched molecular pair (MMP) networks indicated that the fluorination of aromatic rings showed no clear trend toward a positive or negative effect on affinity. Additionally, we propose an in silico workflow (including induced-fit docking, molecular dynamics, quantum polarized ligand docking, and binding free energy calculations based on the Generalized-Born Surface-Area (GBSA) model) to score the fluorine positions in the molecule.

Preclinical characterization of CPL304110 as a potent and selective inhibitor of fibroblast growth factor receptors 1, 2, and 3 for gastric, bladder, and squamous cell lung cancer.

Fibroblast Growth Factor Receptors (FGFRs) are a family of receptor tyrosine kinases expressed on a plethora of cell membranes. They play crucial roles in both embryonic development and adult tissue functions. There is an increasing amount of evidence that FGFR-mediated oncogenesis is mainly related to gene amplification, activating mutations, or translocation in tumors of various histological types. Dysregulation of FGFRs has been implicated in a wide variety of neoplasms, such as bladder, gastric, and lung cancers. Given their functional significance, FGFRs emerge as promising targets for cancer therapy. Here, we introduce CPL304100, an innovative and highly potent FGFR1-3 kinase inhibitor demonstrating excellent in vitro biological activity. Comprehensive analyses encompassed kinase assays, cell line evaluations, PK/PD studies surface plasmon resonance studies, molecular docking, and in vivo testing in mouse xenografts. CPL304110 exhibited a distinctive binding profile to FGFR1/2/3 kinase domains, accompanied by a good safety profile and favorable ADMET parameters. Selective inhibition of tumor cell lines featuring active FGFR signaling was observed, distinguishing it from cell lines lacking FGFR aberrations (FGFR1, 2, and 3). CPL304110 demonstrated efficacy in both FGFR-dependent cell lines and patient-derived tumor xenograft (PDTX) in vivo models. Comparative analyses with FDA-approved FGFR inhibitors, erdafitinib and pemigatinib, revealed certain advantages of CPL304110 in both in vitro and in vivo assessments. Encouraging preclinical results led the way for the initiation of a Phase I clinical trial (01FGFR2018; NCT04149691) to further evaluate CPL304110 as a novel anticancer therapy.

Design, Synthesis, and Development of Pyrazolo[1,5-a]pyrimidine Derivatives as a Novel Series of Selective PI3Kδ Inhibitors: Part II-Benzimidazole Derivatives.

Phosphoinositide 3-kinase (PI3K) is the family of lipid kinases participating in vital cellular processes such as cell proliferation, growth, migration, or cytokines production. Due to the high expression of these proteins in many human cells and their involvement in metabolism regulation, normal embryogenesis, or maintaining glucose homeostasis, the inhibition of PI3K (especially the first class which contains four subunits: α, ß, γ, δ) is considered to be a promising therapeutic strategy for the treatment of inflammatory and autoimmune diseases such as systemic lupus erythematosus (SLE) or multiple sclerosis. In this work, we synthesized a library of benzimidazole derivatives of pyrazolo[1,5-a]pyrimidine representing a collection of new, potent, active, and selective inhibitors of PI3Kδ, displaying IC50 values ranging from 1.892 to 0.018 µM. Among all compounds obtained, CPL302415 (6) showed the highest activity (IC50 value of 18 nM for PI3Kδ), good selectivity (for PI3Kδ relative to other PI3K isoforms: PI3Kα/δ = 79; PI3Kß/δ = 1415; PI3Kγ/δ = 939), and promising physicochemical properties. As a lead compound synthesized on a relatively large scale, this structure is considered a potential future candidate for clinical trials in SLE treatment.

Overcoming undesirable hERG affinity by incorporating fluorine atoms: A case of MAO-B inhibitors derived from 1 H-pyrrolo-[3,2-c]quinolines.

The incorporation of the fluorine motif is a strategy widely applied in drug design for modulating the activity, physicochemical parameters, and metabolic stability of chemical compounds. In this study, we attempted to reduce the affinity for ether-à-go-go-related gene (hERG) channel by introducing fluorine atoms in a group of 1H-pyrrolo[3,2-c]quinolines that are capable of inhibiting monoamine oxidase type B (MAO-B). A series of structural modifications guided by in vitro evaluation of MAO-B inhibition and antitargeting for hERG channels were performed, which led to the identification of 1-(3-chlorobenzyl)-4-(4,4-difluoropiperidin-1-yl)-1H-pyrrolo[3,2-c]quinoline (26). Compound 26 acted as a reversible MAO-B inhibitor exhibiting selectivity over 45 targets, enzymes, transporters, and ion channels, and showed potent glioprotective properties in cultured astrocytes. In addition, the compound demonstrated good metabolic stability in rat liver microsomes assay, a favorable safety profile, and brain permeability. It also displayed procognitive effects in the novel object recognition test in rats and antidepressant-like activity in forced swim test in mice. The findings of the study suggest that reversible MAO-B inhibitors can have potential therapeutic applications in Alzheimer's disease.

Hydrogen Bonds with Fluorine in Ligand-Protein Complexes-the PDB Analysis and Energy Calculations.

Fluorine is a common substituent in medicinal chemistry and is found in up to 50% of the most profitable drugs. In this study, a statistical analysis of the nature, geometry, and frequency of hydrogen bonds (HBs) formed between the aromatic and aliphatic C-F groups of small molecules and biological targets found in the Protein Data Bank (PDB) repository was presented. Interaction energies were calculated for those complexes using three different approaches. The obtained results indicated that the interaction energy of F-containing HBs is determined by the donor-acceptor distance and not by the angles. Moreover, no significant relationship between the energies of HBs with fluorine and the donor type was found, implying that fluorine is a weak HB acceptor for all types of HB donors. However, the statistical analysis of the PDB repository revealed that the most populated geometric parameters of HBs did not match the calculated energetic optima. In a nutshell, HBs containing fluorine are forced to form due to the stronger ligand-receptor neighboring interactions, which make fluorine the "donor's last resort".

Data-Driven Analysis of Fluorination of Ligands of Aminergic G Protein Coupled Receptors.

Currently, G protein-coupled receptors are the targets with the highest number of drugs in many therapeutic areas. Fluorination has become a common strategy in designing highly active biological compounds, as evidenced by the steadily increasing number of newly approved fluorine-containing drugs. Herein, we identified in the ChEMBL database and analysed 1554 target-based FSAR sets (non-fluorinated compounds and their fluorinated analogues) comprising 966 unique non-fluorinated and 2457 unique fluorinated compounds active against 33 different aminergic GPCRs. Although a relatively small number of activity cliffs (defined as a pair of structurally similar compounds showing significant differences of activity -ΔpPot > 1.7) was found in FSAR sets, it is clear that appropriately introduced fluorine can increase ligand potency more than 50-fold. The analysis of matched molecular pairs (MMPs) networks indicated that the fluorination of the aromatic ring showed no clear trend towards a positive or negative effect on affinity; however, a favourable site for a positive potency effect of fluorination was the ortho position. Fluorination of aliphatic fragments more often led to a decrease in biological activity. The results may constitute the rules of thumb for fluorination of aminergic receptor ligands and provide insights into the role of fluorine substitutions in medicinal chemistry.

Influence of Fluorine Substitution on Nonbonding Interactions in Selected Para-Halogeno Anilines.

A series of 4-halogeno aniline derivatives was studied employing combined theoretical and experimental methods (i. e. crystal structure analysis and vibrational spectroscopies). This simplified model system was selected to shed light on the impact of fluorine substitution on the formation of noncovalent interactions such as halogen bonds (XBs) and hydrogen bonds (HBs), which are key interactions in fluorinated/halogenated drug-protein complex formation. Comparative analysis of three previously reported and five newly determined crystal structures indicated that, in most cases, 2-fluoro and 2,6-difluoro substitution of 4-X anilines increases the ability of adjacent amine to form strong N-H⋅⋅⋅N HBs. Additionally, fluorine substituents in the difluorinated derivatives are competitive and attractive HB and XB acceptors and increase the probability of halogen-halogen contacts. A peculiar observation was made for 4-iodoaniline and 2,6-difluoro-4-iodoaniline, which form distinct interaction patterns compared to the corresponding 4-Cl and 4-Br analogs. The observed intramolecular N-H⋅⋅⋅F interactions lead to additional NH bands in the FT-IR spectra.

Tuning the activity of known drugs via the introduction of halogen atoms, a case study of SERT ligands - Fluoxetine and fluvoxamine.

The selective serotonin reuptake inhibitors (SSRIs), acting at the serotonin transporter (SERT), are one of the most widely prescribed antidepressant medications. All five approved SSRIs possess either fluorine or chlorine atoms, and there is a limited number of reports describing their analogs with heavier halogens, i.e., bromine and iodine. To elucidate the role of halogen atoms in the binding of SSRIs to SERT, we designed a series of 22 fluoxetine and fluvoxamine analogs substituted with fluorine, chlorine, bromine, and iodine atoms, differently arranged on the phenyl ring. The obtained biological activity data, supported by a thorough in silico binding mode analysis, allowed the identification of two partners for halogen bond interactions: the backbone carbonyl oxygen atoms of E493 and T497. Additionally, compounds with heavier halogen atoms were found to bind with the SERT via a distinctly different binding mode, a result not presented elsewhere. The subsequent analysis of the prepared XSAR sets showed that E493 and T497 participated in the largest number of formed halogen bonds. The XSAR library analysis led to the synthesis of two of the most active compounds (3,4-diCl-fluoxetine 42, SERT Ki = 5 nM and 3,4-diCl-fluvoxamine 46, SERT Ki = 9 nM, fluoxetine SERT Ki = 31 nM, fluvoxamine SERT Ki = 458 nM). We present an example of the successful use of a rational methodology to analyze binding and design more active compounds by halogen atom introduction. 'XSAR library analysis', a new tool in medicinal chemistry, was instrumental in identifying optimal halogen atom substitution.

How can fluorine directly and indirectly affect the hydrogen bonding in molecular systems? - A case study for monofluoroanilines.

Hydrogen bonds (HBs) directly engaging fluorine has been extensively studied, but the indirect effect of fluorine on adjacent donors and acceptors is poorly understood and still difficult to predict. The indirect and direct effect of the fluorination of aniline on HB patterns observed in monofluoroanilines was studied via experimental (vibrational spectroscopy and crystal structure analysis) and theoretical (ab initio molecular dynamics and electrostatic surface potential) methods. It was found that a fluorine substituent decreases the strength and frequency of N-H⋯N HBs and, at the same time, increases the acidity of CH protons, enhancing the competitiveness of weaker interactions. Additionally, the position of fluorine in the aromatic ring strongly affects the C-F bond length, and a direct intramolecular N-H⋯F HB causes an increase in the N-H bond stability. We also provide a methodology to identify and separate individual HBs concerning the type of donor or acceptor from the ab initio molecular dynamics trajectories.

Imidazopyridine-Based 5-HT6 Receptor Neutral Antagonists: Impact of N1-Benzyl and N1-Phenylsulfonyl Fragments on Different Receptor Conformational States.

G-protein coupled receptors (GPCRs) exist in an equilibrium of multiple conformational states, including different active states, which depend on the nature of the bound ligand. In consequence, different conformational states can initiate specific signal transduction pathways. The study identified compound 7e, which acts as a potent 5-hydroxytryptamine type 6 receptor (5-HT6R) neutral antagonist at Gs and does not impact neurite growth (process controlled by Cdk5). MD simulations highlighted receptor conformational changes for 7e and inverse agonist PZ-1444. In cell-based assays, neutral antagonists of the 5-HT6R (7e and CPPQ), but not inverse agonists (SB-258585, intepirdine, PZ-1444), displayed glioprotective properties against 6-hydroxydopamine-induced and doxorubicin-induced cytotoxicity. These suggest that targeting the activated conformational state of the 5-HT6R with neutral antagonists implicates the protecting properties of astrocytes. Additionally, 7e prevented scopolamine-induced learning deficits in the novel object recognition test in rats. We propose 7e as a probe for further understanding of the functional outcomes of different states of the 5-HT6R.

Rationally designed N-phenylsulfonylindoles as a tool for the analysis of the non-basic 5-HT6R ligands binding mode.

Among all of the monoaminergic receptors, the 5-HT6R has the highest number of non-basic ligands (approximately 5% of compounds stored in 25th version of ChEMBL database have the strongest basic pKa below 5, calculated using the Instant JChem calculator plugin). These compounds, when devoid of a basic nitrogen, exhibit high affinity and remarkable selectivity. Despite a decade of research, no clues have been given for explanation of such an intriguing phenomenon. Here, a series of analogs of four known 5-HT6R ligands, has been rationally designed to approach this issue. For each of the synthesized 42 compounds, a binding affinity for 5-HT6R has been measured, together with a selectivity profile against 5-HT1AR, 5-HT2AR, 5-HT7R and D2R. Performed induced fit docking and molecular dynamics experiments revealed that no particular interaction was responsible for the activity of non-basic compounds. In fact, a plain N-phenylsulfonylindole (1e) was found to possess a moderate (5-HT6R, Ki = 159 nM) affinity. No other monoaminergic receptor has as simple and selective ligand as this one. Thus, it is stated that it binds to the receptor solely based on its conformation and as such, possesses a minimum amount of features, required for binding. Also, any functional group able to form an additional interaction with the receptor increase the binding affinity, like in the case of two highly active non-basic compounds 3e and 5g (5-HT6R, Ki = 65 nM and 38 nM, respectively).

2-Aminoimidazole-based antagonists of the 5-HT6 receptor - A new concept in aminergic GPCR ligand design.

A new strategy in the design of aminergic GPCR ligands is proposed - the use of aromatic, heterocyclic basic moieties in place of the evergreen piperazine or alicyclic and aliphatic amines. This hypothesis has been tested using a benchmark series of 5-HT6R antagonists obtained by coupling variously substituted 2-aminoimidazole moieties to the well established 1-benzenesulfonyl-1H-indoles, which served as the ligands cores. The crystallographic studies revealed that upon protonation, the 2-aminoimidazole fragment triggers a resonance driven conformational change leading to a form of higher affinity. This molecular switch may be responsible for the observed differences in 5-HT6R activity of the studied chemotypes with different amine-like fragments. Considering the multiple functionalization sites of the embedded guanidine fragment, diverse libraries were constructed, and the relationships between the structure and activity, metabolic stability, and solubility were established. Compounds from the N-(1H-imidazol-2-yl)acylamide chemotype (10a-z) exhibited high affinity for 5-HT6R and very high selectivity over 5-HT1A, 5-HT2A, 5-HT7 and D2 receptors (negligible binding), which was attributed to their very weak basicity. The lead compound in the series 4-methyl-5-[1-(naphthalene-1-sulfonyl)-1H-indol-3-yl]-1H-imidazol-2-amine (9i) was shown to reverse the cognitive impairment caused by the administration of scopolamine in rats indicating procognitive potential.

Dual 5-HT6 and D3 Receptor Antagonists in a Group of 1H-Pyrrolo[3,2-c]quinolines with Neuroprotective and Procognitive Activity.

In light of the multifactorial origin of neurodegenerative disorders and some body of evidence indicating that pharmacological blockade of serotonin 5-HT6 and dopamine D3 receptors might be beneficial for cognitive decline, we envisioned (S)-1-[(3-chlorophenyl)sulfonyl]-4-(pyrrolidine-3-yl-amino)-1H-pyrrolo[3,2-c]quinoline (CPPQ), a neutral antagonist of 5-HT6R, as a chemical template for designing dual antagonists of 5-HT6/D3 receptors. As shown by in vitro experiments, supported by quantum chemical calculations and molecular dynamic simulations, introducing alkyl substituents at the pyrrolidine nitrogen of CPPQ, fulfilled structural requirements for simultaneous modulation of 5-HT6 and D3 receptors. The study identified compound 19 ((S)-1-((3-chlorophenyl)sulfonyl)-N-(1-isobutylpyrrolidin-3-yl)-1H-pyrrolo[3,2-c]quinolin-4-amine), which was classified as a dual 5-HT6/D3R antagonist (Ki(5-HT6) = 27 nM, Ki(D3) = 7 nM). Compound 19 behaved as a neutral antagonist at Gs signaling and had no influence on receptor-operated, cyclin-dependent kinase 5 (Cdk5)-dependent neurite growth. In contrast to the well characterized 5-HT6R antagonist intepirdine, compound 19 displayed neuroprotective properties against astrocyte damage induced by doxorubicin, as shown using 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium (MTT) staining to assess cell metabolic activity and lactate dehydrogenase (LDH) release as an index of cell membrane disruption. This feature is of particular importance considering the involvement of loss of homeostatic function of glial cells in the progress of neurodegeneration. Biological results obtained for 19 in in vitro tests, translated into procognitive properties in phencyclidine (PCP)-induced memory decline in the novel object recognition (NOR) task in rats.